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Abstract 

For simple models of the nuclear charge operator, measurements of the 
Coulomb sum and the charge form factor of a nucleus directly determine 
the proton-proton correlations. We examine experimental results obtained 
for few-body nuclei at Bates and Saclay using models of the charge operator 
that include both one- and two-body terms. Previous analyses using one-body 
terms only have failed to reproduce the experimental results. However, we 
find that the same operators which have been used to successfully describe the 
charge form factors also produce substantial agreement with measurements 
of the Coulomb sum. 
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It has long been known that the integrated strength of the longitudinal response function 
measured in inclusive electron scattering (the Coulomb sum rule) is related to the Fourier 
transform of the proton-proton distribution function (PPDF) in the nuclear ground state 
A crucial assumption in obtaining this relation is that the nuclear charge distribution 
arises solely from the protons. As the PPDF is sensitive to the short-range proton-proton 
correlations, its experimental determination can provide direct information on both the 
strength of the correlations in the nuclear medium and the size of the repulsive core in the 
nucleon-nucleon interaction. 

Beck has recently obtained an experimental PPDF from the Bates @ and Saclay Q 
longitudinal data on 3 He. His analysis has shown that a large discrepancy exists between the 
experimental PPDF and that calculated [Q from an essentially exact Faddeev wave function 
H corresponding to a realistic Hamiltonian with the Argonne Vu two-nucleon and Urbana 
VII || three-nucleon interaction models. Specifically, he found that the experimental PPDF 
has a zero at lower momentum transfer and a far greater magnitude in the region of the 
second maximum than the calculated PPDF. It is important to point out that Faddeev 
calculations based on different realistic two-nucleon interactions all give very similar results, 
as reported by Doyle et al. ||. Beck's analysis implies that the experimental PPDF is 
smaller at short distances than the calculated PPDF, thus suggesting that the proton-proton 
interaction has a stronger repulsion than present models would indicate. 

In this letter we analyze the longitudinal response data on 3 H, 3 He, and 4 He obtained at 
Bates [|3|,|10| and on 3 He and 4 He obtained at Saclay [PUTT]. The Coulomb sum is defined as 

where k and u are the momentum and energy transfers, -Rl(/c, uS) is the longitudinal response, 



Ge, p is the proton electric form factor (the Holder parameterization [E| is used in the present 
work), and io e \ is the energy of the recoiling A-nucleon system with Z protons. It can be 
expressed as 

s L (k) = i < o|pi(k) PL (k)|o > -i| < o| Pi (k)|o > | 2 
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where |0 > is the ground state of the nucleus, Fi(k) is the charge form factor normalized 
as = 0) = 1, and a longitudinal-longitudinal distribution function (LLDF) has been 

defined as 
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In this work we assume that the nuclear charge operator pi,(k) consists of one- and two-body 
parts 



p L (k) = p L ,i( k ) +PL,2(k) . 



(4) 



The one-body part includes, in addition to the dominant proton contribution, the neu- 
tron contribution and the Darwin-Foldy and spin-orbit relativistic corrections to the single- 
nucleon charge operator 
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where k 2 = k 2 — (k 2 /2m) 2 is the four-momentum transfer corresponding to the quasielastic 
peak, and Ge,u, Gm,ti an d Gm, p are the neutron electric, neutron magnetic and proton 
magnetic form factors, respectively, evaluated at k 2 . The Darwin-Foldy correction is taken 
into account by the factor 1/(1 + k 2 /4m 2 ) ^ as suggested by Friar |Tj|. Note that because of 
the definition of Sl in Eq. (1), the charge operator in Eq. (4) is divided by Ge, p - 

The two-body part contains contributions associated with pion, p- and tu-meson ex- 
changes, and the pir'f and ujtt^ mechanisms [ TA ] . In the momentum transfer range of interest 
(k ~ 600MeV/c) the pion term is by far the most important. For example, the contributions 



to the A = 3 and 4 charge form factors of the vector meson terms are at least one order of 
magnitude smaller. The pion term is given by 

3i 
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where m n and f n are the pion mass and the 7rNN coupling constant, respectively, with 
/^/47r = .081. The form factor at the 7rNN vertex is chosen to be large (A,r = 2 GeV), as 
suggested by an analysis of the pseudoscalar component of the Argonne vu interaction |L5| . 



F{ and F^ are the Dirac isoscalar and isovector nucleon form factors. The charge operator 
given in Eq. (4) gives an excellent description of the charge form factors of 3 H, 3 He, and 
4 He in calculations based on essentially exact Faddeev (A = 3) |16| and Green's function 
Monte Carlo (GFMC) (A = 4) ]17] wave functions obtained from the Hamiltonian containing 
the Argonne vxa and Urbana VIII interactions. (This Hamiltonian correctly reproduces the 
experimental binding energies of A = 3 and 4 nuclei in Faddeev and GFMC calculations.) 

In order to experimentally determine the LLDF in Eq. (2) it is necessary to measure 
both the charge form factor F^ and the Coulomb sum Sl- For the charge form factors 
of 3 H, 3 He, and 4 He we have used accurate fits to the world data provided to us by Sick 
18fl . As the longitudinal response can be measured only up to some u; max < k by inclusive 



electron scattering, it is necessary to estimate the contribution of the unobserved strength 
for uj > cj max in order to obtain the Coulomb sum. We have assumed that for uj > u max the 
longitudinal response can be parameterized as 

R L (k, u > w max ) = R L (k, u max ; exp) f • (H) 



This form has been suggested by a study of the high w-behavior of the deuteron longitudinal 



response, which can be accurately calculated |T!|. It has been found that for the Argonne 



V\4 interaction, the power a(k) in the deuteron is in the range 3.0-3.5 for k between 200-600 
MeV/c. In the A = 3 and 4 nuclei it is determined by requiring that the energy- weighted 
sum rule W^k) , 



reproduces that calculated as 



W L (k) = ^< 0\p[(k)[H,p L (k)}\0 > ~cu d \ < 0|p L (k)|0 > | 2 , (13) 

by exact Monte Carlo methods. Here H is the Hamiltonian with the Argonne Vu and Urbana 
VIII interactions, and pi is the operator given in Eq. (4). The parameter a is typically found 
to be in the range 2.8-3.5 (2.6-3.1) for the A = 3 (A = 4) nuclei and k = 200-600 MeV/c, 
and does not depend significantly on the value c^Wx chosen. The present analysis differs 
from that reported in Refs. [20,21] in two respects. First, in Ref. [20] the tail contribution 
to SL(k) is estimated by parameterizing the high-cj tail of the response as a sum of two 
decreasing exponentials required to join the data smoothly and to satisfy the calculated 
energy- and energy-square-weighted sum rules. It should be noted that the values for a 
reported above suggest that the energy-square- weighted sum rule may not exist. Secondly, 
the energy-weighted sum rule has been calculated here with a charge operator that includes 
both one- and two-body components rather than the proton contribution only as in Refs. 
[20,21]. The two-body components (predominantly that associated with pion exchange) 
lead to an enhancement of W L {k) of 10% (6.0%), 8.6% (4.3%), 7.5% (3.4%) in 3 H ( 3 He) and 
9.1%, 7.8% and 7.4% in 4 He at k = 300, 400, 500 MeV/c, respectively. The dominant kinetic 
energy contribution, which is exactly given by k 2 /2m when only protons are included in p^, 
is little affected by the relativistic corrections and two-body terms. However, the latter 
enhance the leading interaction contributions associated with isospin-exchange spin and 
tensor components. 



As a consequence of these differences, the present analysis yields values for S^k) that 
are slightly larger (~2%) than those published in Ref. [20] for 3 H and 3 He, the only data 
for which a comparison is possible. The final analysis of the Bates data on 4 He, published 
in Ref. [10], has given a separated longitudinal response that is somewhat smaller in the 
quasielastic peak than that used in Ref. [20] to obtain S^k). 

The experimental LLDF obtained for 3 H, 3 He, and 4 He are compared with theory in 
Figs. 1-3. The errors in the experimental LLDF are dominated by those in the Coulomb 
sum. The latter has two sources: the first, from the measured portion of S^k), denoted 
as Si{k; exp); the second, from the tail contribution, denoted as ^(fcjtail). The error on 
Si(k] exp) has been estimated by adding in quadrature the random errors on the measured 
longitudinal response function and by further assuming the systematic error to be as large 
as the random error so obtained. The error on Sl (k; tail) has been estimated by assuming 
it to be given by Sl(/c; tail) x AR^k, u) mauX )/ Ri(k, u max ), where ARl is the experimental 
error on R L (k,u) at to = cu max (typically ~20-30% of R^k, tu max ))- 

The theoretical curves in Figs. 1-3 have been obtained by exact Monte Carlo evaluation 
of the expectation value in Eq. (3). We have used exact Faddeev (A=3) and GFMC (A=4) 
wave functions, again corresponding to the Argonne v 14 plus Urbana VIII interaction models. 
The LLDF obtained from the Bates and Saclay data on 3 He and 4 He are consistent with 
each other, within errors, and are in good agreement with the results of calculations in 
which both one- and two-body terms are included in the charge operator. In particular, 
the position of the zero and magnitude of the second maximum are well reproduced by 
these calculations. The results obtained by neglecting the contributions due to the two- 
body terms or by keeping only the proton contributions in pi are in poor agreement with 
the data: the zero is shifted to higher fc's and the strength in the second maximum is 
greatly underestimated, as found by Beck |2|. We also note that in the calculation of 3 H, 
3 He and 4 He charge form factors, inclusion of the two-body components in pi is crucial 
for correctly reproducing the experimental data in the diffraction minimum region. The 
minima are located at k ~ 630 MeV/c in the 3 He and 4 He charge form factors. However, 
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they occur at significantly lower momentum, k ~ 370 MeV/c, in the 3 He and 4 He LLDF, 
thus enhancing the importance of the relativistic and meson-exchange corrections at low 
momentum transfers. 

In 3 H the LLDF calculated in the approximation in which only protons are considered 
vanishes identically. However, the experimental LLDF extracted from the Bates data is 
different from that obtained in the full calculation. This discrepancy is also found in the 
Coulomb sum rule: the experimental Si(k) (including the tail contribution) is larger by ~ 
10% in the k = 350-500 MeV/c range than the theoretical one. However, the 3 H experimental 
charge form factor is well reproduced by the present theory. 

To summarize, the LLDF has been calculated in the A = 3 and 4 nuclei with exact 
Faddeev and GFMC wave functions obtained from a realistic Hamiltonian containing the 
Argonne v u two-nucleon and Urbana VIII three-nucleon interaction models. The charge 
operator has been taken to include, in addition to the dominant proton contribution, also 
the neutron contribution, the Darwin-Foldy and spin-orbit relativistic corrections, and two- 
body terms associated with meson exchanges. Within this framework good agreement has 
been obtained between the calculated and experimental 3 He and 4 He LLDF. However, large 
discrepancies remain between the calculated and experimental 3 H LLDF. 

The present 3 He and 4 He results indicate that, because of the complicated nature of the 
coupling between a longitudinal virtual photon and the nucleus (even at low momentum 
transfers), the LLDF extracted from inclusive electron scattering data cannot provide di- 
rect information on the strength of the proton-proton repulsive interaction at short range. 
Therefore, the experimental evidence based on inclusive longitudinal data for proton-proton 
short-range correlations remains elusive. 
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FIGURES 

FIG. 1. Experimental and theoretical longitudinal-longitudinal distribution functions in 3 He. 
Circles (squares) denote Bates (Saclay) data; solid symbols denote negative values. The curves 
labelled proton, 1-body, and full show theoretical results obtained from the Faddeev wave function 
by including in pL the proton, one-body, and one- plus two-body contributions, respectively. 

FIG. 2. Same as in Fig.l but for 4 He with theoretical results from the GFMC wave function. 

FIG. 3. Same as in Fig.l but for 3 H. 
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